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Abstract 

We propose a method to determine the phase motion of a complex amplitude 
in three body heavy meson decays. We show that the phase variation of a 
complex amplitude can be directly revealed through the interference observed 
in the Dalitz Plot region where it crosses with a well established resonant state. 
This method could be applied to the decays D + —> 7r _ 7r + 7r + and D + —> 
K~tt + tt + , to determine whether the low mass states a and ac, suggested by 
E791, have phase motions compatible with resonances. 
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Recently the Fermilab experiment E791 with Dalitz plot analysis showed strong 
evidence for the existence of two light and broad scalar resonances. One resonance 
is compatible with the isoscalar meson cr, and is observed in Cabbibo suppressed 
decay D + —> 7T~'n + 7i + [if]. The other, its strange partner, the meson k seen in the 
D + —» K~ti + ti + decay 0. These results opened a new experimental window for 
understanding scalar meson spectroscopy 0. 

In these analyses, each possible resonance amplitude is represented by Breit- 
Wigner functions multiplied by angular distributions associated with the spin of the 
resonance. The various contributions are combined in a coherent sum with complex 
coefficients that are extracted from maximum likelihood fits to the data. The abso¬ 
lute value of the coefficients are related to the relative fraction of each contribution 
and the phases takes into account final state interaction (FSI) between the resonance 
and the bachelor pion HI. To fit E791 data, it was necessary to include extra, not 
previously observed scalar resonances. For the new states, modeled by Breit Wigner 
amplitudes, they measure mass and width. The D + —> 7r _ 7r + 7r + decay required a 
scalar with mass 478^3 ±17 MeV/c 1 2 and width 324^ ±21 MeV/c 2 . The high statis¬ 
tics sample of D + —> K~h + ti + needed a scalar with mass of 797 ± 19 ± 43 MeV/c 2 
and width 410 ± 43 ± 87 MeV/c 2 for a good confidence level fit. In both cases, the 
extra contribution is dominant, accounting for approximately half of the decay. Due 
the importance of these mesons in many areas of particle and nuclear physics, it is 
desirable to be able to have confirmation with a direct observation of their phase 
motion, without having to assume a priori the Breit-Wigner phase variation |dj. 

Partial wave analyses (PWA) have been used successfully in hadron-hadron scat¬ 
tering to observe the phase motion of complex amplitudes. For example, the LASS 
experiment 0 did a partial wave analysis of the K~n + spectrum in K~p —» K~n + n 
interaction with transverse momentum less than 0.2 (GeV/c) 2 . The low momentum 
transfer assures first that one 7r + meson exchange between the incoming K~ and 
the proton is dominant, and second that there is no important contribution to final 
state interactions between the K~n + system and the outgoing neutron. This is a 
physical situation completely different from charm three body decay, where all phase 
space is taken into account, including low and high momentum transfer between the 
resonance state and the bachelor particle. The strength of the resonance-bachelor 
particle interaction in charm decay is given by the FSI phase. In fact, experimental 
results show a tendency of large values for the strong phase, thus non negligible FSI 
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The method proposed here regards the Dalitz plot just in the region where a well 
established resonance, represented by a Breit-Wigner, for example, in S 12 spectrum @, 

1 This phase is considered constant because it depends only on the total energy of the system, i.e. 
the heavy meson mass. 

2 In the three body decay D —> ijk, Sij and Sjk are two of the three possible Dalitz plot variables, 
that is the square of the invariant mass of the particles i plus j and j plus k respectively. 
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crosses the complex amplitude under study in the Si 3 variable. In Figure [I] we show 
a Monte Carlo simulation of a Dalitz plot of two scalar resonances, one the / 0 (980) 
in the variable S 12 and the other a broad and low mass resonance compatible with 
the a in S 13 . We show that the interference created in this particular region of the 
phase space, can reveal the unknown phase variation of the amplitude as a function of 
s 13 . The method can be applied to retrieve slow variation similar to the one observed 
in Kit elastic scattering at low invariant mass |J, as well as rapidly moving ones 
expected for resonant amplitudes. 



Figure 1: Dalitz plot distribution for two crossing scalar mesons, /o(980) in the vari¬ 
able S 12 and a in S 13 . The left and right hashed regions correspond to | A(m^— e, S 13 ) | 2 
and to | M(? 7 Tq + e, S 13 ) | 2 respectively 

For simplicity we suppose that the only contributions we have in the Dalitz plot 
are the amplitude we want to study in s 13 and a scalar state represented by a constant 
width Breit-Wigner in Si 2 . This is not a limitation of the method, it being straight¬ 
forward to include known effects like extra states populating the region as well as 
angular momentum distributions and varying width Breit-Wigner functions. In this 
simple case, the total amplitude is given by: 

A(s 12 , S13) = a R BW(s 12 ) + a s sinS(si 3 ) e*^ Sl3)+7) ( 1 ) 

7 represents the constant FSI phase difference between the two body scattering am¬ 
plitude and the bachelor particle. a R contains all constant or almost constant pa¬ 
rameters, like form factors and the magnitude. a s e lS ^ Sl3 ^sind(s 13 ) is the most general 
form to represent a complex amplitude, a s can be considered constant in the range 
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of the invariant mass region of interest. The relativistic Breit-Wigner function used 
in many Dalitz plot analyses is given by: 




m 0 T 0 

«12 - ml + im 0 T 0 


( 2 ) 


mo To are mass and the width of the resonance. The square modulus of the above 
amplitude can be written as: 


I - 40 i 2 ,si 3 ) | 2 = a R | BW{s 12 ) I 2 +a 2 sm 2 h(si 3 ) + 

2a R a s m 0 Tosinti (* 13 ) ^^ _ m 2 ) C 0 S ( 5 ( Sl3 ) + 7 ) + m 0 r 0 sm(h(si 3 ) + 7 ) ( 3 ) 

O 12 •m5) 2 + mgT5 

I 11 this case, | £TfV(si 2 ) | J is a symmetric function and thus £TfV(si 2 = ml + e) 
equals BW(s' u = ml — e), where e is a small Si 2 interval. Consequently the difference 
of the amplitudes square takes the simple form: 

| A(ml + e, si 3 ) | 2 - | A{m 2 0 - e, s 13 ) | 2 = ^^—^-sinS(s 13 )cos(5(s 13 ) + 7 ). (4) 

e + m 0 l 0 

| ^4(mg — e, S 13 ) | 2 and | A(ml + e, S 13 ) | 2 are taken from data as becomes clear from 
Figure |T|. We can rewrite Eq. 4 as: 

A | A | 2 =| -4(mo + e, si 3 ) | 2 - | A(ml - e, si 3 ) | 2 = C(sin(28(si 3 ) + 7 ) - sin'y). (5) 

C is equal to 4a s aijemoT 0 /(e 2 + mgTg). 7 being a constant, A | A | 2 directly reflects 
the behavior of J(si 3 ). A constant A | A | 2 would imply constant h(si 3 ), this would 
be the case of non-resonant contribution. The same way a slow phase motion will 
produce a slowly varying A | A | 2 and a full resonance phase motion produces a clear 
signature in A | ^4. | 2 with the presence of zero, maxim and minim values. Next we 
discuss how to extract the constant phase 7 and <5(si 3 ) for two situations. We start 
with the case of a rapidly moving phase like a resonance where we extract direct 
information on h(si 3 ). Then we deal with slow phase variation, for which we use as 
the example the LASS Kit low mass phase behavior. 

To study a resonant amplitude in the A | A | 2 distribution, we produced a fast 
simulation of 10K events with a hypothetical decay D + —> 7 r“ 7 r + 7 r + with only a broad 
and low mass scalar resonance cr in S 13 and the well known resonance /o(980) in Si 2 
and strong phase difference between the two of 2.8 rad, Figure [I] shows the relevant 
portion of the Dalitz plot. To generate the / 0 (980) we used the constant width 
Breit-Wigner (Eq. 2). For the low mass wide cr, the Breit-Wigner parameterization 
produces undesirable threshold features. We could, in principle, introduce form- 
factors to avoid this problem, but this would bring unnecessary complexity to the 
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example. Instead ,we parameterized the low mass resonance as sinS(si 3 )e lS ^ Sl3 \ like 
the amplitude we used in Eq. 1. with £(si 3 ) = 7 r / (1+ exp (a(si 3 — ml)), the 
parameter a is associated with the a width. We used a = —15 GeV~ 2 . With this 
choice, we have a complete n phase variation around the a mass, as would behave a 
Breit-Wigner far from the threshold, and the mass plot given by | sin5(si 3 )e* <5 ^ Sl3 ' ) | 2 
reproduces the overall characteristics of the a resonance. 




Figure 2 : a) A j A | 2 distributions for resonance a like interfering with /o(980) gener¬ 
ated with a 10 K fast-MC events. b)h(si 3 ) distribution got from A | A | 2 distribution 
presented in a). The dashed line is the function used as input to the fast-MC. 

The A | A j J distribution for the fast-MC produced with the above parameters 
is presented in Figure Qa. We can see a strong variation of this distribution with 
the presence of three zeros: one at the threshold implicit from Eq. 5, the second at 
s 33 ~ 0.3 associated with 5 = 37 t /2 — 7 and the third at s i3 ~ 0.55 corresponding 
to S = 7 r. Thus, it is possible to observe in a simple way, a clear signature of the 
presence of large phase variation associated with the number of zeros in A | A | 2 . 

We can extract directly the phase motion h(si 3 ) noticing that the maximum and 
minimum values of A | A | 2 corresponds respectively to the sin(28(si 3 ) + 7 ) = 1 and 
sm(25(s 13 ) + 7 ) = —1 conditions. These and the A | A | 2 value, contents of bins 3 
and 6 of the plot Figure [2]a provide us two equations that can be used to get 7 and 
C. Then we invert Equation 5 and get from each value of A | A \ 2 one value of h(si 3 ). 
Figure |2|b compares h(si 3 ) obtained in this way with the functional form (dashed) 
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Figure 3: 7 distribution for 1000 fast-MC experiments, phase 2.8 rad. 


used to generate the fast-MC sample. To evaluate error in J(si 3 ) (Figure ||b), and the 
ability of this procedure to measure the phase 7 we repeated the fast-MC experiment 
1000 times. The result is shown in Figure [|. We are able to recover the input value 
( 2.8 rad) and the rms of the distribution is consistent with the individual experiment 
errors. 

We have presented an ideal situation in our example, with constant width, and 
isolated scalar resonances. In a realistic case, the width dependence introduces an 
asymmetry in S 13 , the same occurs with respect to resonant amplitudes with spin 
1 or 2 and even form factors [TO. These effects can be included in Equation 3 in 


a complete way, since we know very well the functional form of a Breit-Wigner, 
the width dependence with s 12 , and the spin amplitude. However, the simplified 
amplitude for A | A | 2 still works if we are willing to sacrifice the measurement 
of 7 . In fact, we observed that many of these realistic effects are absorbed by the 
parameter 7 , keeping the <5 (si 3 ) distribution almost unchanged if we include or not 
the realistic effects. Since with this method, one is primarily interested on the direct 
measurement of the phase motion, and since the usual Dalitz plot method has been the 
best technique to get the FSI phase, one can use A | A | 2 without much elaboration 
in most cases. 

The other case that we wish to discuss is when the complex amplitude is not 
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a resonance, but has only a small phase variation. A good example is the scalar 
smooth phase variation at low mass obtained by the LASS 0 collaboration in Kir 
scattering. The LASS phase variation is given by cot5 = 1/d.q + e.q/2 ||, where q is 
the momentum in the K~ti + center of mass, d = 4.03±1.72 GeV -1 and e = 1.29±0.89 
GeV -1 are fit parameters. We included this amplitude as a function of si 3 in a 
hypothetical three body charm decay D + K~7i + n + with the well known resonance 
Aq( 1430) in the crossing channel, s i2 . Figure £|a shows the A | A | 2 behavior, given 
by Eq. 5, for several values of 7 between 0° to 300°. While all curves show slow 
variation of A | A | 2 , we notice the clear distinction among the 7 values. 

In this case we can not extract directly the phase motion since we have not maxima 
and minima, as we have for the resonance example. Instead we have to use a functional 
form to fit data and get the shape. The fit function should be monotonic with a 
constant asymptotic behavior. We used for the 5 (si 3 ) form the same function we 
used to generated the resonance in the preceding example, that is <f(si 3 ) = ph / ( 1 + 
exp (a(si 3 — t)), where ph, a and t are free fit parameters. With this simple function 
we are able to retrieve the <5 (si 3 ) behavior and 7 values. In Figure f|b we show the 
LASS phase behavior [5J, input to the exercise and the dashed line the functional 
form we get with our fit function. 
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Figure 4: (a)A | A | 2 distributions for different 7 values, (b) LASS low mass phase 
motion used as input (dotted line) and obtained fitting A j A | 2 distribution using 
the functional form of <5 (si 3 ) given in text (dashed line). 

We have presented here a method to observe phase motion of a complex amplitude 
in a Dalitz plot of three body decays. We showed that the difference of amplitudes, 
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A | A | 2 , around the central mass squared of a well known resonance is a good way 
to observe this variation. We also showed that one can determine rapid or even slow 
phase variation with this method. We used as example two scalar resonances, the a 
and the /o(980) in D + —> 7r - 7r + 7r + decay. This is the simplest possible situation, 
however one could use also the crossing between the a and the /2(1270) resonance 
including the well known angular distribution. The advantage of using / 2 ( 1270) in 
real data is the large statistics in the D + —> 7r - 7r + 7r + decay, associated with the size 
of the branching fraction and the effective area of the crossing region. Nevertheless, 
we believe that both / 0 (980) and / 2 ( 1270) could be used to measure the phase motion 
of the low 7r + 7r _ invariant mass region. The analysis of the k resonance in D + —> 
K~tt + tt + decay, should be slightly more complicated because of the proximity of the 
vector resonance K* (890) to the new scalar state. However the K* (890) would have 
to be incorporated in Eq. 1. 

We thank Professors Jeff Appel, Hans Giienter Dosch, Wolfgang Ochs and Alberto 
Reis for suggestions and important comments. 
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